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Abstract. We show, by means of a perturbative weakly nonlinear anal- 
ysis, that the axisymmetric magneto-rotational instability (MRI) in 
a magnetic Taylor-Couette (mTC) flow in a thin-gap gives rise, for 
very small magnetic Prandtl numbers (Pm), to a real Ginzburg-Landau 
equation for the disturbance amplitude. The saturation amplitude As 
is found to scale in this regime as 

VL, with 1/2 < 5 < 2/3 (depend- 
ing on the boundary conditions adopted). The asymptotic results are 
shown to comply with numerical calculations performed by using a 
spectral code. They suggest that the transport due to the nonlinearly 
developed MRI may be vanishingly small for Pm ^ 1. 



1 Introduction 

• 1.1 General 



, Asymptotic approaches to nonlinear stellar pulsation were pioneered by Robert 
' Buchler and his collaborators (see this volume). Amplitude equations, which 
• • , capture the essential dynamics near the instability threshold, result from such 
' approaches and may usually be derived using singular perturbation theory (see 
, Marie- Jo Goupil, this volume). In this contribution we report on the results of 
applying a perturbative asymptotic analysis to an instability that has acquired, 
5^ , in the past 15 years or so, a paramount importance in the quest to understand 
angular momentum transport in accretion disks, namely the magneto-rotational 
instability (MRI). Even though the system studied here is rather different from 
pulsating stars, the ideas and techniques are quite similar. One important differ- 
ence is that while the amplitude equations for stellar pulsation are sets of ordinary 
differential equations, here we must take recourse to partial differential equation(s) 
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(PDE), in which the amplitude (usually called in this case envelope) is a function 
of time (as in the ordinary amplitude equations) and space. The application of 
singular perturbation theory requires thus the use of slow variation in time as well 
as in space. In this way spatio-temporal slow dynamics (i.e., patterns) is captured 
by a generic PDE - the real Ginzburg-Landau Equation, rGLE, in our case. 

The rudiments of singular perturbation theory are well explained in the book 
by Bender & Orszag (|1999p . while a detailed account on problems dealing with 
the dynamics on the center manifold can be found in Manneville (|1990p . It seems 
that these powerful analytical and semi-analytical techniques have not yet found 
their way to enough astrophysical applications. In the early 1980s Robert Buch- 
ler started an ambitious program of introducing such a new approach to stellar 
pulsation. I was fortunate to be his postdoc then, but have moved on to other 
topics. Robert and his postdocs (most of which contributed to this volume) have 
pursued this program and during this conference one could learn just how much 
has been achieved in the understanding of stellar pulsation, following asymptotic 
approaches, combined with numerical and phenomenological ones. It is quite un- 
fortunate that "driving forces" of Robert's magnitude have unfortunately been 
quite rare in other branches of astrophysics. 

1.2 The MRI - background 

The linear MRI has been known for almost 50 years (Velikhov ll9591 Chandrasekhar 
I1960p : Rayleigh stable rotating Couette flows of conducting fluids are destabilized 
in the presence of a vertical magnetic field if dfl^/dr < (angular velocity de- 
creasing outward). However, the MRI only acquired importance to astrophysics 
after the influential work of Balbus & Hawley (|1991|) . who demonstrated its vi- 
ability to cylindrical accretion disks (for any sufficiently weak field). The linear 
analysis was later supplemented by nonlinear numerical simulations (albeit of a 
small, shearing-box (SB), segment of the disk). Enhanced transport (conceivably 
turbulent) is necessary for accretion to proceed in disks, found in a variety of as- 
trophysical settings - from protostars to active galactic nuclei. As is well known, 
Keplerian rotation law is hydrodynamically linearly stable and thus the MRI has 
been widely accepted as an attractive solution for enhanced transport (see the re- 
views by Balbus & Hawlev [1 9981 and Balbus [?003p . even though some questions on 
the nature of MRI-driven turbulence still remain (see, e.g., Branderburg|2005p, in 
particular in view of some recent numerical studies (Pessah, Chan & Psaltis 2 0071 
Fromang & Papaloizou I2007| Fromang et al. I2007p . Because of the MRFs im- 
portance and some of its outstanding unresolved issues, several groups have quite 
recently undertaken projects to investigate the instability under laboratory con- 
ditions. Additionally, numerical simulations specially designed for experimental 
setups have been conducted. It appears, however, quite impossible to make defi- 
nite deductions from these experiments and simulations (see, e.g., Ji et al. 120011 
Noguchi et al. 120021 Sisan et al. 120041 Stefani et al. 120061 Stnd references therein). 

We shall report here on the result of a weakly nonlinear analysis of the MRI 
near threshold, for a dissipative mTC flow. We have done the analysis for the above 
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simplified (relatively to an accretion disk) problem for two types of boundary con- 
ditions (see Umurhan, Menou & Regev [^071 hereafter UMR and Umurhan, Regev 
& Menou 1^0071 hereafter URM, for a detailed account). This kind of approach 
is important because the viability of the MRI as the driver of turbulence and 
angular-momentum transport relies on understanding its nonlinear development 
and saturation. By complementing numerical simulations with analytical methods 
useful physical insight can be expected. To facilitate an analytical approach we 
made a number of simplifying assumptions so as to make the system amenable 
to well-known asymptotic methods (see, e.g., Cross and Hohenberg 120031 Regev 
I2006p for the derivation of a nonlinear envelope equation, valid near the linear 
instability threshold. 



2 Linear theory 

2.1 Equations 

The hydromagnetic equations in cylindrical coordinates (Chandrasckhar ll96ip are 
applied to the neighborhood of a representative radial point (ro) in a mTC setup 
with an imposed constant background vertical magnetic field. The steady base 
flow has only a constant vertical magnetic field, B — Bqz, and a velocity of 
the form V = U{x)y. In this base state the velocity has a linear shear profile 
U{x) = —qflox, representing an azimuthal fiow about a point tq, that rotates with 
a rate fto, defined from the differential rotation law il{r) cx f2o(f'/?'o)~' {q = 3/2 for 
Keplerian rotation) . The total pressure in the base state (divided by the constant 

density), 11 = po^^ (^P + , is a constant and thus its gradient is zero. 

This base flow is disturbed by 3-D perturbations on the magnetic fleld b = 
{bx, by,bz), as well as on the velocity - u = (u^, Uy,Uz), and on the total pressure - 
w. We consider only axisymmetric disturbances, i.e. perturbations with structure 
only in the x and z directions. This results, after non-dimensionalization, in the 
following set of equations, given here in the rotating frame: 

^ - 217oz X u - gf^o^i^y - Cb • Vb - CBo<9^b = -Vtu + ^^-V^u, (2.1) 
^ - b • Vu + gf^ofo.y - Soa.u = :^V2b, (2.2) 

V • u EE d^Ux + d^Uz =0, V • b = d^b,., + d.b^ = 0. (2.3) 

The Cartesian coordinates x,y, z correspond to the radial (shear- wise) , azimuthal 
(stream-wise) and vertical directions, respectively, and since axisymmetry is as- 
sumed V = yidx + idz and = 9^ + 9f . Lengths have been non-dimensionalized 
by L (« the gap size), time by the local rotation rate f2o (tildes denote dimensional 
quantities). Because the dimensional rotation rate of the box (about the central 
object) is rJo = ^^oz, the non-dimensional quantity is simply equal to 1, but we 
keep it to flag the Coriolis terms. Velocities have been scaled by J^qL and the mag- 
netic field by its background value Bq. Thus = 1 as well, but again, we leave it 
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in the equation set for later convenience. The hydrodynamic pressure is scaled by 
PqL^CIq and the magnetic one by 3^/(8^). The non-dimensional perturbation tu 
of the total pressure divided by the density (which is equal to 1 in non-dimensional 
units), which survives the spatial derivatives, is thus given hy ru = p + Ci|bp, 
where p is the hydrodynamic pressure perturbation. The non-dimensional param- 
eter C = B'q/{A-kpdUqL'^) = V^/V"^ is the Cowling number, measuring the relative 
importance of the magnetic pressure to the hydrodynamical one. It is equal to 
the inverse square of the typical Alfven number {Va is the typical Alfven speed). 
The Cowling number appears in the non-linear equations, together with the two 
Reynolds numbers: TZ = Q.qL'^ /v and TZm ^ ^oL'^/v, where i) and fj are, respec- 
tively, the microscopic viscosity and magnetic resistivity of the fluid. We shall also 
see that the magnetic Prandtl number, given as Vm = T^m/T^, plays an important 
role in the nonlinear evolution of this system. 

It is useful to rewrite the above equations of motion in terms of more convenient 
dependent variables. Because the flow is incompressible and y-independent, the 
radial and vertical velocities can be expressed in terms of the streamfunction, 
that is, {ux,Uz) = {dz'^ ,—dx'^). Also, since the magnetic field is source free, 
one can similarly express its vertical and radial components in terms of the flux 
function, <3>, that is, (bx,bz) = (9^$, -dx'^) (see URM for details). The system is 
supplement by appropriate boundary condition on the thin cylindrical gap " walls" , 
as detailed in UMR and URM and we also point out that the stress relevant for 
angular momentum transport, S, is 

S = Sr -|- Em, Sr = UxUy, Em = —Cbxby. 

which consists of Sr and Em, the Reynolds (hydrodynamic) and Maxwell stresses, 
expressing the velocity and magnetic field disturbance correlations, respectively. 



2.2 Linear stability 

Linearization of (2.1-2.3), for perturbations of the form cx g^t+^kxx+ik^z ^ gives rise 
to the dispersion relation, 

V{s; kx,kz,Vm,S,C,q) = aos'^ + ais^ + a28^ + a^s -h 04 = 0, 

where all the a, are functions of kx, kz and the other parameters, but we shall 
write explicitly only the coefficient that will be used in what follows, 



C 

a4 = ^ 



klC{CklVm + klS^f + K^SW^kl - 2qS% , (2.4) 



where the notation = 2(2 — q), fc^ = k'^ + k^ is introduced. 

For given values of the parameters (call them pj) there will be four distinct 
modes. The linear theory in various limits for this problem has been discussed in 
numerous publications and we shall not elaborate upon it any further. Rather, 
we focus on situations where the most unstable mode (of the four) is marginal 
(critical) for some 3jij a given value of k^ = K. We fix K in order to focus on 
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Fig. 1. Summary of linear theory. This example is for C = 0.08, Vm ~ 0.001, q = 3/2, 
and the fundamental mode, (a) Growth rates, Re{s), as a function of wavenumber k for 
three values of TZm- (b) Solid line depicts those values of TZm and k where Re{s) — 0. 
The shaded region shows unstable modes. The locations of fc = kcrit = Q and TZm = 
7?,m(crit) = Rm are shown. 



the marginal vertical dynamics within the thin gap (actually a rotating channel, 
see below). Marginality (s = 0) implies the vanishing of the real coefficient a4{kz) 
(as expressed above) and its derivative with respect to kz at some kz — Q: 

ai{kz=Q;K,pj)^Q, ^{kz = Q; K,pj) = 0. (2.5) 

The second condition and (|2.4p yield 

q[q^S^{CQ^ - 4q) + [K^ + 02)2g2c52(52 + k") + 
{K^ + Q'^fCS'^i&CVn.Q'' + k') + 

2{K^ + Q'^fVmC'^S'^ + 5{K^ + QY'Pm^C^] = 0- (2-6) 

This equation, together with a4{kz — Q,K:pi) = 0, can be solved for S and Q. 
The general expressions for S{K, Vm,C, q) and Q{K, Vm,C, q) are lengthy but their 
asymptotic forms, to O (Pm) (for Va\ ^ 1), are simple. 



^16 Cq{2~q)K .,2 '^q-CK^ 
^= {2q-CK^) ' ^="^2^70^ ^'-^^ 

If CK^ > 2q the solutions arc not physically meaningful, while CK^ = 2q corre- 
sponds to the ideal MRI limit. 

Our system is a rotating thin channel, whose walls are at x = 0,Tr/K. All 
quantities are vertically periodic on a scale Lz commensurate with integer multiples 
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of 2tt/Q. As long as ^/Q, the limit of a vertically extended system (and 

thus a continuous spectrum of vertical modes) may be effected. Linear theory, as 
discussed above, is summarized in Figure [1] 

3 Weakly nonlinear analysis 

We follow the fluid into instability by tuning the vertical background field away 
from the steady state, i.e. Bq ^ 1 — e^A, where e <C 1 and A is an O (1) control 
parameter. This means that we are now in a position to apply procedures of 
singular perturbation theory to this problem, by employing the multiple-scale (in 
z and t) method (e.g, Bender & Orszag [1999p . It facilitates, by imposing suitable 
solvability conditions at each expansion order of the calculation in order to prevent 
a breakdown in the solutions, a derivation of an envelope equation for the unstable 
mode. Thus for any dependent fluid quantity F{x, z, t) we assume 



The fact that the x and z components of the velocity and magnetic field pertur- 
bations can be derived from a streamfunction, ^{x, z) and magnetic flux function 
$(a;, z) reduces the number of relevant dependent variables F to four [uy and by 
are the additional two). These four dependent variables will also be used in the 
spectral numerical calculation (see below). 

For the lowest e order of the equations, resulting from substituting the expan- 
sions into the original PDE and collecting same order terms, we make the Ansatz 
Fi(x,z,t) = FiA(e2;, e^f)e''^^ sinifa; -I- c.c, where Fi is a constant (according to 
the variable in question), and where the envelope function A (an arbitrary con- 
stant amplitude in linear theory) is now allowed to have weak space (on scale 
Z = ez) and time (on scale T = e^t) dependencies. Because this system is tenth 
order in x-derivatives, a sufficient number of conditions must be specified at the 
edges, and the Ansatz has to obey them. In UMR we chose mathematically expe- 
dient boundary conditions, which allowed fully analytical treatment for the limit 
'Pm ^ 1 and here we shall report, in some detail, only on the results of that work. 
We have also tried some different boundary conditions (see URM) and in that case 
the rGLE coefficients had to be calculated numerically. The qualitative behavior 
of the saturation amplitude was found to be quite similar (see below). 

3.1 Real Ginzubrg-Landau equation 

The end result of the asymptotic procedure procedure is the well-known real 
Ginzburg-Landau Equation (rGLE) which, for Vm ^ 1, is 



F{x, z, t) = eFi{x, z, t) + e^F2{x, z,t) + ■ ■ ■ . 



(3.1) 



dxA = XA 



1 



\A\A^ + DdlA. 



(3.2) 



Vn.C 



We shall consider only the magnitude (real part) of A. In a ID gradient system 
like this one phase dynamics can only give rise to wavelength modulations and 
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consequently our main results are in no way affected by it. The real envelope is 
A = y/^A, A = (X, and the coefficients, for q = 3/2 are 



3 5S^ - 185^ - 32 + 2{S^ + WWS^ + 1 
4' 5(52 + l)(4V52TT-3) 



(3.3) 



D = 6^ ^ ^ 3.4 

53(4V5^n-3) ^ ' 

^ 3SVS^ - 65 

452 + 1 + X/52TT' ■ 

where 5(i4r) is as given in (|2.7p . For 5 ^ f , (i.e. as one approaches the ideal MRI 
limit), these simplify to ^ = 15/16, ( ~ 3/4, D — 3/2; and in general they remain 
O (1) quantities for all reasonable values of 5. 

3.2 Scaling of angular momentum transport 

Contributions to the total angular momentum transport (J = Jh + Jb) due to the 
hydrodynamic (Jh) and magnetic correlations (Jb) are, to leading order, 



The envelope A is found by solving the rGLE, which is a well-studied system 
(see, e.g., Regev [20061 for a summary and references). It has 3 steady uniform 
solutions in ID (here, the vertical): an unstable state ^ = 0, and 2 stable states 
A = ±As, where A^ = |C(5)|7'inC (note that C(5) is an O (1) quantity). As is also 
the saturation amplitude, because the system will develop towards it. 

Setting A ^ As in (13. 6p and in the expression for Ey (the total disturbance en- 
ergy, see UMR) , followed by some algebra, reveals that the angular momentum flux 
in the saturated state is, to leading order in Vm and e, J = e'^\({S)\'^{S)VmCS^^ , 
while in the expression for the energy one term is independent of Vmi Ev = 
e'^C^C'^{S)l32iS) + 0{V,^). The key resuffs of this analysis (see URM) are the 
following scalings (with the magnetic Prandtl number) 



V^m^Ev-i^o, J-^m (3.7) 

for Vm 1) where Eq is a constant, independent of T'm 0- For fixed resistivity 
this imphes J '-^ TZ~^ . 



^For the boundary conditions used in URM we numerically got As oc 
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Time (in box orbits) 



Fig. 2. Ev (panel a) and J (panel b) as a function of time from numerical calculation. 
The different lines are labeled by the value of the Reynolds number TZ The diamonds 
in panel (b) show the scaling predicted by our asymptotic analysis, predicting also a 
constant final value of the disturbance energy, as apparent in panel (a). One "box orbit" 
= 2ti in nondimensional units. 



We have also performed fully numerical calculations, using a 2-D spectral code 
to solve the original nonlinear equations, in the streamfunction and magnetic flux 
function formulation, near MRI threshold. The code implemented a Fourier- 
Galerkin expansion in 64 x 64 modes in each of the four independent physical 
variables, i.e. F = ^ Fnm(i) siniir„a;e*'^"^+c.c., where F = {ip , Uy , ^ , by)'^ and 
where Fnni(i) is the time-dependent amplitude of the particular Fourier-Galerkin 
mode in question (denoted by the indices n,m). We typically started with white 
noise initial conditions on Fnin(O) at a level of 0.1 the energy of the background 
shear. Because of space limitations we display only a representative plot of runs 
made with S = 5.0, C = 0.05, (i.e. 7^,„ fixed), q = 3/2 and = 0.2 for a few 
successively increasing values of TZ (Figure [2]). With this value of e the fastest 
growing linear mode has a growth rate of ~ 0.065, and thus fully developed ideal 
MRI can not be expected. Ey saturates at a constant independent of TZ (for small 
enough T'm), while J saturates at values that scale as TZ~^. Thus the asymptotic 
analysis fits very well the fully numerical results after the system saturates. 



4 Summary 

The numerical and asymptotic solutions developed here show that for T'^ <C 1, 
it is the azimuthal velocity perturbation, arising from the O (e^) term, that be- 
comes dominant in the saturated state. It appears to be the primary agent in 
the nonlinear saturation of the MRI in the channel, acting anisotropically so as to 
modify the shear profile and results in a non-diagonal stress component (relevant 
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for angular momentum transport). Our analysis is complementary to the study 
of Knobloch and Julien (|2005p . who have performed an asymptotic MRI analysis, 
but for a developed state, far from marginality. 

The trends predicted by this simplified model are not qualitatively altered by 
different boundary conditions, the value of S changing only by a small amount. 
Qualitatively similar scalings, but with somewhat different values of 5 > 0, have 
very recently been reported by Lesur & Longaretti (|2007p and Fromang et al. 
P007p . however in those SB simulations was not taken to be vanishingly small. 
It should be remarked that conventional accretion disks do have ■pm ^ 1, but the 
numerical SB simulations cannot faithfully treat very small such Prandtl numbers, 
for numerical reasons. Further analytical investigations of the kind reported here 
should contribute to better physical understanding of the MRI saturation. . 
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